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ABSTRACT 

We examine whether the gravitational shear-intrinsic cllipticity (GI) correlation function of the 
luminous red galaxies (LRGs) can be modeled with the distribution function of a misalignment angle 
advocated recently by Okumura et al. For this purpose, we have accurately measured the GI cor- 
relation for the LRGs in the Data Release 6 (DR6) of the Sloan Digital Sky Survey (SDSS), which 
confirms the results of Hirata et al. who used the DR4 data. By comparing the GI correlation func- 
tions in the simulation and in the observation, we find that the GI correlation can be modeled in the 
current ACDM model if the misalignment follows a Gaussian distribution with a zero mean and a 
typical misalignment angle a$ = 34.9+2 1 degrees. We also find a correlation between the axis ratios 
and intrinsic alignments of LRGs. This effect should be taken into account in theoretical modeling of 
the GI and intrinsic ellipticity-ellipticity correlations for weak lensing surveys. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: halos — gravitational 
lensing — large-scale structure of universe — methods: statistical 



1. INTRODUCTION 

Weak gravitational lensing by large-scale structure is 
a promising tool to directly probe matter distribution 
in the universe. The main source of contaminations for 
weak lensing observations comes from two types of in- 
trinsic alignments: the ellipticity correlation of source 
galaxies with each other [intrinsic ellipticity-ellipticity 
(II) correlation] and the ellipticity correlation of lense 
galaxies with the surrounding matter distribution [grav- 
itational shear-intrinsic ellipticity (GI) correlation]. In- 
vestigating these intrinsic alignments is also important 
for galaxy formation studies. 

There has been much theoret i cal work (e.g.. 
Heavens et all 120001: |Croft fc Metzlerl 120001 : iLee k Pen! 
2000t ICatelan et all 
serva tional work 
2001 



20011: Uind 120021) a s well as ob 



iHirata et all l2004t 



Pen et alJ l2000t iBrown et all 



Mand clbaum et al.1 [2006; 



IHevmans et al l 12004 
lOkumura et al.l l2009f) 

which attempted to estimate the II correlation. In 
lOkumura et al.l (|2009l ). we have determined the II cor- 
relation of luminous red galaxies (LRGs) accurately to 
a large scale (~ 30 h^ 1 Mpc) using the Dat a Release 6 
of th e Sloan Digital Sky Survey (SDSS DR6; lYork et all 
2000). From this measurement, we further gave a 
constraint on the misalignment angle of typically 35° 
(26° on average) between the central giant elliptical 
galaxies and their host dark matter halos using the 
halo occu pation distributio n (H OD) approach (e.g., 
I Jing et all 119981: IZheng et al.ll2008f) . This misalignment 
is found to be in agreement with the misalignment 
between the inner dark matter and the host halos by 
iFaltenbacher et al.l ((2009) using an Y-body simulation, 
though it is unknown if this is a coincidence because 
the light distribution of LRGs is not necessarily to be 
the same as the inner dark matter distribution of halos. 
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Nevertheless, the observational constraint does provide 
very useful clues to understanding the contaminations 
on weak lensing su rveys as well as to stu dying formation 
of giant ellipticals ((Okumura "et~aT1l2009h . 

Given the distribution of the misalignment angles, we 
can also predict the GI of LRGs with the surrounding 
dark matter distribution for current cosmological models 
by combining an Y-body simulation and the HOD mod- 
eling. Such predictions can be co mpared with the ob - 
ser vations of the GI such as t hose of lHirata et al.1 (|2007f ) 
andlFaltcnbac her et alJ (|2009t ) . The comparisons will fur- 
ther test if the distribution of the misalignment angle is 
correct, and will furthermore tell us how to model the 
GI contaminations generally for future weak lensing sur- 
veys. This approach is complementary to recently pro- 
posed statistical approaches which attempt to eliminate 
or minimize the GI effect in weak lensing surveys pure! 
based on observational samp les (Hirata & Scliakl 1200 
[Kml [20051: IHevmans et~aT| l200aTBridle fc Kin? 120071: 
Joachimi fc Schneider! 12008: Zhang 20081: iHui & Zhand 
120081) . We will present such a study in this Letter. 

2. SDSS LUMINOUS RED GALAXY SAMPLE 

Ou r HOD model p arameters used in Section [4] are 
from iSeo et al.l (l2008h which is based on the fitting of 
IZheng et al.l (|2008|) to the proje cted correlation function 
of LRGs of IZehavi et all 0001). Whe n we started the 
work, we tri e d to use the observational data of GI of 
IHirata et alJ {2007), but later we found that the pro- 
jected correlation function presented in their paper is 
slightly smaller (about 20 % in low-z and 15 % in high- 
z subsamples) than that in Ze havi et al.l (|2005h . While 
we do not know if the samples they used are slightly 
different, i t is a pparent that the HOD parameters of 
Seo et all (120081) c annot be applied to the sample of 
Hirata et alJ (2007). Therefore, we decided to measure 



the GI of LRGs using the SPSS DR6 data (lYork et all 
20001 lEisenstein et al.1 120011 : lAdelman-McCarthv et"aTI 
2008( 1. The LR G sample is the same as that used in our 



previous paper (|Okumura et al.ll2009h . There are 83,773 
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Fig. 1. — Projected galaxy-galaxy correlation function for the 
LRG sample at 0, 16 < z < 0.36. For co mparison, the data of 
IZehavi etall (120051) and lHirata et al.l fl2007f ) are also plotted. 

LRGs in the redshift range 0.16 < z < 0.47. 

Here we need to model the radial and angular selec- 
tion functions of the observational sample. We build 
the radial selection function by a spline fit with a Gaus- 
sian smoothing of the redshift distribution of observed 
LRGs. The angular selection function is constructed 
using the angular survey mas k provided by the Value 
Added Galaxy Catalog (VAGG lBlanton et We 
assumed that our sample has the same survey geometry 
as that of the VAGC and excluded all the LRGs which 
are not overlapped. Then the number of LRGs is reduced 
to 78,758. Finally we id entify 73,935 cen t ral LR Gs using 
the criteria proposed by iReid fc Spergell (|2008| ) . 

We measure the projected correlation function from 
the LRG sample, w gg (r p ) = J (, gg (r p , TI)cflT, where 
£, gg (r p ,H) is the galaxy autocorrelation function as a 
function of separations perpendicular (r p ) and paral- 
lel (IT) to the line of sight and is measured using the 
lLandy fc Szalavl (|1993| ) estimator. Figure Q] shows the 
co mparison of our mea sured w gg with the previous work 
by IZehavi et al . ( 2005). Very good agreement between 
the two studies confirmed that our selection functions 
work well on the scales of interests in this work. 

To measure the GI function of LRGs, we also need 
information on their shapes. The ellipticity of galax- 
ies is determined by the SDSS photometric pipeline 
called photo and defined as the ellipticity of the 25 ma g 
arcsec^ 2 isophote in the r band (Stought on et al.| [2002). 
The point-spread f unction has been co rrected when mea- 
suring the shapes (jLupton et al.ll200l[) . The components 
of the ellipticity are defined as 



1 - q 2 f cos 2/3 
sin 2/3 



1 



(1) 



where q is the ratio of minor and major axes and /3 is the 
position angle of the ellipticity from the north celestial 
pole to east. 

3. GRAVITATIONAL SHEAR-INTRINSIC ELLIPTICITY 
CORRELATION OF LRGS 

To estimat e the GI correlation, we ad opt t he formalism 
develo ped bvlMandclba um et al. ( 2006) an dlHirata et al.1 
(|2007fh The generalized lLandv fc Szalavl dl993j) estima- 
tor is used for estimating the GI correlation function, 



S+(D-R) 
RR 



(2) 



where RR is the normalized counts of random-random 
pairs in a particular bin in the space of (r p ,H). S+D is 
the sum over all pairs of the + component of shear in 

(r PI n), 

e +(j'l«) 



S+D 
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where the ellipticity component of jih LRG, e + (j\i), is 
redefined relative to the direction to the zth LRG and 
thus corresponds to the elongation along the direction, 
and 1 Z = 1 — cr 2 N = 1 ~ ( e +) is t ne shear responsivity 
(e.g.. iBernstein fe Jarvisl l2002f ) and K « 0.947 for our 
LRG sample. S + R is calculated likewise using a cata- 
log of randomly distributed points in the survey region. 
Finally the projected GI correlation function is obtained 
by doing the projection along the radial direction, 



w. 



g+ 



(r P ) 



-n„ 



£„ + (r p ,n)dII. 



(4) 



We adopt n max = 80 h~ x Mpc, but changing n max from 
60 to 100 h^ 1 Mpc does not significantly change w g+ (r p ) 
for r p < 100 h~ x Mpc. Positive w g+ (r p ) means that 
the major axes of LRGs tend to point toward over- 
densities at a transverse scale r p . w g + is related to 
the density-intrinsic ellipticity correlation function w$+ 
throug h the galaxy bias b Q as w a + = b a ws-\- on large 
scales (jMandelbaum et al.lf200l IHirata et al.ll2007l see 



also Section [5]). 



w 



9X 



(r p ) is also calculated in the same 



way and used for a test of systematics because w gx (r p ) 
should be zero on all scales. 

The resulting GI correlation functions for the observed 
central LRGs are shown in Figure [U The error bars 
represent la errors estimated from 93 jackknifed realiza- 
tions. Here the axis ratio q in equation {T]) is set to be 
zero. This is equivalent to the assumption that a galaxy 
is a line along its major axis. In this case the shear 
responsivity becomes 72. = 1 — (cos 2 2/3} = 0.5. Clear 
detection of the GI correlation w g + can be seen up to 
large scales. Note th at when our result is compared with 
the previous work of IHirata et a l. (2007), the difference 
of the q factor and thus of 1Z must be taken into account. 
The prescription of q = and its relationship with that 
of q > will be discussed in Section 14.31 From the plots 



of w gx in Figure [2] we confirm that all the values are 
consistent with zero at the 2a confidence level. We thus 
discuss the GI correlation only in terms of w g + in the 
following analysis. 

4. COMPARISON TO MODEL PREDICTIONS 

4.1. Modeled Gravitational Shear-Intrinsic Ellipticity 
Correlation Function 

To make model predictions for the GI correlation, 
we fol low the same methodology as in lOkumura et al.l 
((20091) . We use a halo catalog constructed from a high- 
resolution cosmological simulation wit h 1024 3 particle s 
in a cubic box of side 1200 hr 1 Mpc ping et all 120071) . 
Central galaxies are assigned to the simulated halos us- 
ing the best-fit HOD p arameters for LR Gs found by 
iSeo et all (|2008l )(see also lZheng et al.ll2008l ). The result- 
ing fraction of mock central LRGs is 93.7% and we use 
only the centrals in order to compare with our observa- 
tion. 
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Fig. 2. — Projected GI correlation functions from the SDSS 
central LRGs, w g +(r p ) and w gx (r p ). The dashed and solid lines 
are respectively the model predictions for to s + obtained from the 
mock LRGs with erg = 0° and 35°. Mixed logarithmic and linear 
scalings are used for the vertical axis. The horizontal axis at the 
top shows the corresponding angular scale when the transverse 
separations are located at z = 0.5. 

We consider halos to have triaxial shapes (|Jing fc Sutol 
2002) . The two components of the ellipticity of each halo 
are estimated from the second moments of the projected 
mass distribution (e.g- ICroft fc Metzlerll2000l) 

V e ' x ) J-xx + Lyy V LLx y ) 

where = i ^ N XiXj and N is the number of particles 
in a halo. Then the GI correlation function of halos is 
measured in the same way as that of LRGs, where the 
value of q is assumed to be zero again and thus 1Z = 0.5. 

First we assume that all central galaxies are completely 
aligned with their parent dark matter halos. The GI cor- 
relation function of central galaxies is then calculated and 
shown in Figure [2l In order to refine the statistics, we 
averaged over seven mock samples with different random 
seeds for assigning LRGs to dark halos. Interestingly, the 
GI correlation function of the mock LRGs, when they are 
assumed to be aligned completely with their host halos, 
has the same shape as but is about twice as high as the 
observ ation. Similar results were shown for the II corre- 
lation (jOkumura et al.ll2009t ). 

4.2. Constraints on Misalignment 

In this subsection, we consider the case in which the 
major axis of each central galaxy is misaligned with that 
of its host halo and give a constraint on the misalignment 
angle by comparing the observed GI corr elation function 
w a + w ith its model prediction. Following lOkumura et all 
(2009) we assume that the misalignment angle 9 between 
the major axes of central LRGs and their host halos fol- 
lows a Gaussian function with a zero mean and a width 
ag, where ag is the typical misalignment angle. We arti- 
ficially assign misalignment to the position angle of each 
mock central LRG relative to its host halo according to 
the Gaussian function. For each chosen value of ag and 
each LRG mock sample, we generate nine misaligned 
LRG samples by choosing different random seeds. Our 
model prediction for each ag is thus calculated by aver- 
aging over 7 x 9 = 63 misaligned samples. 



In comparing the observational data with the model 
prediction, v 2 statistics are calculated in the range of 
20° < ag < 50°. In this analysis we use the seven data 
points of w g ^(r p ) shown in Figure [2l while there is one 
free parameter, ag] thus the degree of freedom is 6. The 
covariance matrix estimated using 93 jackknifed subsam- 
ples is used for the calculation of x 2 - 

The fits of the observed GI correlation function w g + 
to the model prediction give a tight constraint on the 
misalignment parameter, ag = 34.9+2? degrees (68% 
confidence level) , which corresponds to a mean misalign- 
ment angle of 26°. 9. This is in very good agreement 
with our previous work on the II correlation which gave 
ag — 35. 4+33 degrees. The constraint from the GI corre- 
lation is tighter than that from the II correlation because 
the GI correlation is better determined. The model pre- 
diction of w g + with ag = 35° is shown in Figure O 

4.3. Correlation of the LRG shape and its orientation 

The misalignment parameter was constrained with the 
assumption of q = 0, i.e., we considered the orientation 
of the LRGs relative to their spatial distribution only. If 
there is no correlation between the shape of the LRGs 
and their orientation, we can use the misalignment an- 
gle distribution to model the GI of LRGs even when 
the shape is inclu ded (i.e., GI in weak lensing studies; 
iHirata et al.ll2007l ). In order to see if this correlation ex- 
ists, we define a normalized GI correlation function w g + 
as 

w g + (r p ; q) = / \ + q q2 \ w 9+ ( r p ; l) > ( 6 ) 

where / \ is the value averaged over all objects in 

the sample (it is 0.29 for observed LRGs and 0.42 for 
mock host halos of LRGs), and w g +(r pi q) is the same 
as equation ([4]) except the q dependence is included. If 
there is no correlation between axis ratios and orienta- 
tions, we expect w g+ {r p ; q = <7lrg) = w g+ (r p ; q = 0) for 
observed LRGs and w g+ (r p ; q = q mock ) = w g+ (r p ; q = 0) 
for halos. Here we neglect the factor of the shear re- 
sponsivity, 1/272., so equation ^ just corresponds to the 



replacement of 1 / 21Z by 




-l 

in equation ([3]). 



In Figure [3] we show w g+ (r p ; 0) of the observed LRGs 
and mock LRGs with ag = 35° , which are the same data 
as those in Figure [2] because w g +(r p ;0) = w g +(r p ;0). 
Next we calculate w g +(r p ; q m ock) for the mock LRGs 
with g m ock determined from the dark matter distribu- 
tion within the host halos, which is shown as the dashed 
line in Figure [3] These results indicate that there ex- 
ists a correlation between the shape and orientation of 
the dark matter halos, and this correlation leads to an 
increase of ~ 15% in the normalized GI correlation. If 
the shapes of the halos were uncorrelated with their ori- 
entations, we expect that this w g +(r p ; q m0 ck) would be 
equal to w g+ (r p ;0). This is indeed confirmed in the fig- 
ure, where the dotted line, the correlation when the val- 
ues Qmock (not the orientation) of host halos are shuffled 
randomly, is completely overlapped with the solid line. 
Almost the same amount of increase is seen for the GIs of 
the observed LRGs, as shown by the squares and the cir- 
cles in the figure, which implies that there is a correlation 
of the shapes and orientation in the observation. 
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Fig. 3. — Normalized GI correlation functions iZ> 9 +(r p ; q). All 
the lines are the model predictions with ag = 35° . The dotted blue 
line is the one with q moc k shuffled randomly. The circles/squares 
have been respectively offset in the horizontally negative/positive 
direction for clarity. 

This correlation should be taken into account when one 
models the GI and II effects in theory by using the distri- 
bution function of the misalignment angle. It should also 
be considered if w g +(r p ;q = (?lrg) is used to constrain 
the distribution function of the misalignment angle. If 
the correlation is neglected, the misalignment angle could 
be slightly underestimated. This is the reason why only 
the or ientation of the LRGs are used in I Okumura et all 
( 2009) and in the current work to constrain the misalign- 
ment angle. 

5. DISCUSSION 

In this Letter, we examined whether the GI correla- 
tion of LRGs can be modeled with the d istribution func- 
tion o f misalignment angle advocated by I Okumura et aT] 
(2009) based on the II correlation. For this purpose, 
we have accurately measured the GI correlation for the 
L RGs in the SP SS DR6, which also confirms the results 
of lHirata et al.l (|2007T ) who used the DR4 data. By com- 
paring the GI correlations in the simulation and in the 
observation, we found that the GI correlation can be 
modeled in the current ACDM model, if the misalign- 
ment follows a Gaussian distribution with a zero mean 
and a typical misalignment angle ag — 34.9^2 ^ degrees. 



This constraint on ag is in excellent agreement with the 
previous work, ag — 35At_f^ degrees, based on the II cor- 
relation. The constraint on ag is tighter in this Letter, 
because the GI correlation is better determined than the 
II correlation in the observation. Furthermore, the good 
agreement of the observed and theoretical GI functions 
further lends nontrivial support for the ACDM scenar- 
ios and for the distribution function of the misalignment 
angle. 

We have found a correlation between the axis ratios 
and intrinsic alignments of LRGs. If the correlation is 
neglected, one would underestimate the GI correlation 
(in the case of q ^ 0) by ~ 15% if the shape q and ori- 
entation 8 of the LRGs are randomly chosen from their 
distribution functions. This effect should be taken into 
account in theoretical modeling of the GI and II correla- 
tions for weak lensing surveys. 

These results have profound implications both for fu- 
ture weak lensing surveys and for studying the formation 
of giant elliptical galaxies. For weak lensing surveys, the 
relevant quantity is the correlation function between the 
mass overdensity a nd the intri nsic ellipticit y , w$+ . As 
iMandelbaum et alJ (|2006f ) and iHirata et all (|2007f ) rea- 
soned, these two quantities are simply related through 
the galaxy bias as w g+ = b g w$+. We can also measure 
w$ + directly in our simulation. Our result supports that 
b g in this simple relation is about 2 for all the scales 
explored here and for our chosen ACDM model. 
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